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Executive  Summary 

Our  long-term  goal  is  the  development  of  automated  tools  for  computing  fluid  flow  in 
complicated  geometries.  In  the  last  decade,  we  have  developed  both  algorithms  and 
software  for  mesh  generation  and  adaptation,  surface  preparation,  multi-grid  coarse  grid 
generation,  and  a  highly  parallel  steady-state  inviscid  flow  solver.  In  combination  with 
other  tools  (in  particular  the  ability  to  generate  high  quality  triangulations  directly  from 
CAD  using  CAPRI,  developed  by  Bob  Haimes  and  Michael  Aftosmis),  complete 
simulations  can  be  performed  with  greatly  reduced  set  up  time.  In  the  last  two  years  a 
significant  advance  is  the  development  of  an  adjoint  module  for  CartSD  by  Aftosmis  and 
Nemec.  This  is  used  to  automatically  drive  the  mesh  adaptation,  and  produce  an  error 
estimate  for  an  objection  function  such  a  lift  or  drag. 

A  specific  goal  of  this  project  was  to  extend  Cartesian  embedded  boundary  methods  to 
solve  the  compressible  viscous  Navier-Stokes  equations.  The  first  step  in  this  direction 
was  to  develop  a  discretization  for  the  viscous  terms.  This  was  done  for  two  dimensional 
problems,  and  several  test  cases  from  the  literature  were  studied.  Future  work  will 
include  extensions  to  three  dimensions,  turbulence  modeling,  and  a  numerical  method  or 
model  for  handling  the  boundary  layer  without  resorting  to  the  prohibitively  expensive 
use  of  isotropic  refinement.  This  work  is  in  collaboration  with  Michael  Aftosmis  and 
Marian  Nemec  of  NASA  Ames. 

Summary  of  Research 

The  CartSD  flow  solver  was  extended  to  include  viscous  terms  in  two  dimensions.  The 
first  step  in  modifying  the  steady-state  inviscid  code  to  simulate  the  Navier-Stokes 
equations  was  to  develop  discretizations  for  the  second  derivative  terms  at  cut  cells  and  at 
mesh  interfaces.  We  examined  four  different  stencils  for  use  at  cut  cells  in  an  elliptic 
model  problem.  We  chose  the  best  one  to  implement  in  the  Navier  Stokes  code.  We 
found  it  necessary  to  use  a  quadratic  reconstruction  in  the  normal  direction  at  the  cut 
cells,  but  we  retained  the  linear  gradient  for  tangential  reconstruction.  We  also 
experimented  with  using  a  second  order  re-centering  at  the  mesh  interfaces  in  the 
direction  normal  to  the  interface.  Doing  this  only  in  the  one  direction  was  not  sufficient 
to  reduce  the  overall  error,  and  we  are  sticking  with  the  linear  reconstruction  in  these 
parts  of  the  mesh.  This  development  work  was  in  two  space  dimensions. 


We  experimented  with  several  test  cases  from  the  literature.  First  we  ran  experiments  for 
a  Blasius  boundary  layer  profile  with  a  non-coordinate  aligned  grid.  As  was  done  by 
Coirier,  we  compared  results  with  a  flat  plate  oriented  with  the  mesh  and  oriented  at  an 
angle  to  the  mesh.  The  Blasius  profiles  are  shown  in  Figure  1,  and  some  results  for  skin 
friction  are  shown  in  Figure  2. 


Figure  1  Boundary  layer  solution  in  self-similar  variables  with  plate  at  several 
angles  to  the  grid,  compared  to  Blasius  solution  in  red. 
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Figure  2  Computed  skin  friction  for  the  boundary  layer  compared  to  the  Blasius 
results.  Left,  from  Coirier’s  thesis  [1];  right,  our  results  (using  the  more  common 
log-log  plot.) 


Another  test  case  is  flow  around  a  NACA0012  at  0  angle  of  attack,  Re  5000,  Mach 
number  .5  The  magnitude  of  the  velocity  if  shown  in  Figure  3.  A  plot  of  Cp  and  skin 
friction  is  shown  in  Figure  4.  The  separation  point  computed  was  81.7,  well  within  the 
range  of  values  computed  for  this  test  case. 


Figure  3  Velocity  magnitude  is  shown  around  a  NACA  0012,  with  Re  5000,  zero 
angle  of  attack,  and  Mach  number  .5 


In  figure  4  left  we  show  the  pressure  coefficient,  and  on  the  right  the  skin,  showing  the 
upper  and  lower  surface.  When  the  mesh  is  not  resolved  and  one  can  see  the  results 
staircasing  through  the  Cartesian  grid.  Most  of  this  calculation  is  well  resolved,  with 
only  a  hint  of  a  problem  at  the  leading  edge  (which  is  hard  to  see  unless  the  figures  are 
blown  up). 


Figure  4  Cp  (left)  and  skin  friction  (right)  for  the  results  on  Figure  3  for  flow 
around  a  NACA  0012. 

A  revised  multigrid  acceleration  was  developed  for  viscous  problems.  The  theory 
suggests  that  higher  order  derivatives  need  a  better  prolongation  than  the  piecewise 
constant  approach  we  used  for  inviscid  flow.  We  developed  a  linear  prolongation 
algorithm,  which  also  seems  to  help  with  robustness.  To  be  able  to  evaluate  viscous 
terms  on  the  coarser  grids  they  need  to  have  some  geometry.  For  example,  the  coarse 
grids  need  to  know  the  cell  centroid,  surface  centroid,  cut  face  information,  and  other 
values,  to  be  able  to  compute  gradients  on  the  coarser  grids.  This  involved  a  substantial 
amount  of  work  developing  the  coarse  mesh  generator.  These  terms  were  not  needed  to 
compute  inviscid  flow,  since  the  coarse  grids  use  only  a  first  order  method  in  this  case, 
and  no  coarser  grid  gradients  are  used.  Figure  5  shows  a  convergence  plot  for  the 
NACA0012  viscous  case  above.  A  full  multigrid  scheme  is  used,  starting  on  the  coarsest 
level  and  using  five  grid  levels. 


Figure  5  Multigrid  convergence  for  the  computation  in  Figure  3  for  flow  around  a 
NACA  0012.  This  case  converges  in  about  200  iterations  on  the  flnest  grid. 


More  infrastructure  will  need  to  be  developed.  Our  Cartesian  mesh  generator  already  has 
an  anisotropic  option  [2],  although  it  has  not  been  pushed  very  far  or  even  fully 
debugged.  We  developed  a  new  multigrid  mesh  coarsening  algorithm  that  uses  a  tree  to 
represent  the  mesh  instead  of  space-filling  curves  to  generate  the  coarse  meshes.  This 
was  necessary  since  space-filling  curves  do  not  retain  their  nice  properties  on 
anisotropically-refined  meshes.  The  tree-based  approach  will  also  allow  for  the  option  of 
variable  mesh  coarsening,  for  example  either  semi-coarsening,  or  more  aggressive 
coarsening  when  mesh  interfaces  slow  things  down.  This  work  was  done  with  my 
Masters  student  Greg  Tumolo,  now  graduated.  Most  important  will  be  the  development 
of  a  boundary  layer  model,  to  reduce  the  prohibitive  cost  of  isotropic  refinement  at  the 
boundary.  This  will  be  the  focus  of  future  research. 
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